Spontaneous Phase Separation in a Vibro-fluidized Granular Layer 
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A new type of two phase coexistence in a vibrated granular material is reported in the limit where energy 
transfer between the particles and the bottom wall occurs due to discrete collisions. A horizontally homogeneous 
bed develops inhomogeneities in which a dense and a dilute phase separate and coexist. The dilute region 
resembles a gas, while the dense region is like a liquid and a velocity distribution very different from a Gaussian. 



Granular materials occur in diverse sizes and shapes. This 
leads to a variety of complex patterns in collective motions 
of grains, which are of relevance to predicting natural phe- 
nomena as well as for process design in industry [|l|] . Though 
the interactions between grains are complex in nature, many 
complex patterns are reproduced by simple microscopic in- 
teractions, such as hard sphere inelastic collisions. Of recent 
interest are granular materials subjected to vertical vibrations. 
When the amplitude of bottom wall vibrations is compara- 
ble to the bed height, the top surface becomes unstable and 
forms a periodic wavy pattern with peaks and troughs . An 
unusual type of excitation is the "oscillons" reported in 
where the excitation is localized. In general it is argued that 
[Q] the various patterns observed in a vibrated granular layers 
are a result of interaction between standing waves. Several of 
these patterns — squares, stripes, hexagons, kinks, etc. — have 
been recovered have by a simple model of interacting waves 

Eggers [|6|] has simulated a vibrated bed under gravity with 
two compartments separated by a wall which allows particles 
from one side to the other through a small slit at a given height. 
For vigorous shaking, the populations in the two halves re- 
main equal, but at a lesser rate of shaking, the bed separates 
into stable coexistence of a high and a low density phase. It 
is also shown that such a stationary state can be captured by 
a flux balance of particles from one side to the other calcu- 
lated by the simple kinetic theory of . Occurrence of large 
scale inhomogeneities in industrial vibrated beds is reported 
in [|^, which seem to be similar in nature to a two phase coex- 
istence, but the instability disappears when the surrounding air 
medium is evacuated. Phase separations have been observed 
in two dimensional horizontal plane simulations and ex- 
periments [pQ|]. The particles are restricted to move only in 
the horizontal two dimensional plane, and the coexistence is 
between two phases whose static structure are an ordered (col- 
lapsed) phase and a liquid-like phase (which is referred to as 
"gas-like" in [0). For a mono-layer of particles, the the ratio 
r of the maximum acceleration of the wall and the accelera- 
tion due to gravity in [ |iq[ ] is less than unity. 

A new type of phase coexistence phenomenon observed in 
two dimensional computer simulations of a vertically vibrated 
granular material is reported here. The parameter regime 
where this is observed differs from that of previous studies 
in two respects. 

1 . These are observed in a low density regime, where the 
area fraction of the particles is small, i.e., the mean free 



path is large compared to the particle diameter. This is 
in contrast to previous studies [ p^ where the mean free 
path less than the particle diameter. 

2. In the present case, the frequency of the base vibration 
is large compared to the frequency of oscillation of the 
bed, and amplitude is small compared to the mean free 
path of the particles. Consequently, the amplitude and 
frequency of the base are not relevant parameters in de- 
termining the state of the bed as in Q, and the only 
relevant parameter is the energy transfer from the base 
to the particles in discrete particle collisions. 

Since we are interested in the dilute limit, it is useful to 
characterize the system by the parameters relevant for a dilute 
vibro-fluidized bed [Q] in the limit where the area fraction is 
small and the coefficient of restitution is close to 1. The rel- 
evant parameters are Ncr, which gives the number of mono- 
layers of particles at rest, and the parameter e = {Uq/T), 
which is a measure of the ratio of the dissipation over time 
scales comparable to the time between collisions and the aver- 
age energy of a particle. Here N is the number of particles per 
unit width, a is the particle diameter, Uq is the mean square 
velocity of the vibrating surface and T, the 'granular temper- 
ature', is the mean square velocity of the particles (the mass 
of a particle is set equal to 1). It was shown that in the limit 
(1 — e^) <C 1, the parameter e is given by 
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The area fraction of particles u in the dilute limit is expo- 
nentially decaying in the vertical direction, and is given by 
where = n N g a'^ / (4T) is the area frac- 
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tion of particles at the bottom of the bed. The area fraction 
at the bottom is kept at a constant value = 0.077 in the 
simulations. 

The other parameters usually used for characterizing pat- 
tern formation in vibrated granular systems are not relevant 
for the present simulations. 

1 . The ratio of characteristic frequency of vertical oscilla- 
tion of the bed {g/VT), and frequency co^ of the bot- 
tom wall oscillations, is maintained at a constant value 
of 0.01 in the simulations, indicating that there are no 
correlated collisions of the bed with the bottom wall. 
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2. The ratio of accelerations F = {U^^uOw/g), which is 
proportional to the inverse of the ratio of frequencies, 
varies in a range 10^-10"^. 

3. The ratio of the amplitude a^^ of the bottom wall vibra- 
tions to the mean free path varies in the range 5 x 10~^- 
0.04, so that the wall does not induce any mean motion 
in the bed. 

In this parameter regime, the dynamics of the bed is coupled 
to that of the bottom wall only through the energy transfer 
to the particles (by random collisions with the wall), and the 
mean motion of the bed is not correlated to the dynamics of 
the bottom wall. Therefore, the periodicity of the bottom wall 
vibration does not induce a periodic motion in the granular 
layer. This is confirmed by the observation that the instabili- 
ties are observed even when the bottom wall is kept stationary, 
but the wall velocity at collision is selected randomly with a 
probability distribution identical to that for the vibrating wall 

The simulation cell is a box containing N circular discs 
per unit width vibrated at the bottom. The gravitational ac- 
celeration force g acts vertically downwards. There are two 
side walls separated by a width and the top is open. Inter- 
particle collisions are inelastic, and are characterized by a con- 
stant coefficient of restitution e. For the sake of simplicity, the 
interaction of the particles with the three walls is assumed to 
be elastic. 

In the present study, we consider the effects varying Ncr 
and e; the other parameters are kept constant at the values 
mentioned above. The simulations show that at low values 
of the parameter e, less than about 0.3, there is a fluidized 
state which is uniform in the horizontal direction. As e is 
increased above 0.3 at constant Ncr < 5, there is a sepa- 
ration of coexisting regions of high and low densities. (For 
Ncr > 5, the system forms convection rolls, which will be 
reported elsewhere). This can be seen in Fig. [l|, which shows 
the coexistence of a dense and a dilute region. After the for- 
mation of two coexisting regions, there is no further variation 
in the average densities of the two regions, though there is a 
variation in the position of the interface between the two re- 
gions. Changing the width of the cell does not alter the length 
scale of the density fluctuation, but simply alters the number 
of dense and dilute regions as in Fig. ^, which shows the varia- 
tion of the scaled local density fluctuation N{x) = Nx/N —1 
for two different cell widths. (It can also be shown that there 
is no phase separation for smaller widths, giving one possible 
reason why such instabilities were not observed earlier.) The 
static structure of the two phases can be represented by the 
pair distribution function g{r), which is shown in Fig. ||. Here 
we see that the dilute region is like a gas phase with uncor- 
rected particle positions and the dense phase is like a liquid 
phase showing correlations up to a few particle diameters. 

Another salient feature of this coexistence is the form of the 
velocity distribution function in the dense phase. Fig. ^ shows 
the velocity distribution functions in the dilute and the dense 
regions. A Gaussian curve with the same temperature is also 
plotted along each for reference. It is seen that in the dilute 
region, the distribution function is very close to the Gaussian 



FIG. 1 : Exploded view of snapshots of regions near the left and the 
right wall. A coexistence of the dilute and dense regions can be 
seen near the right wall. The parameters used for this simulation 
are: iVcr = 1, e = 0.75, and z/° = 0.077. The horizontal width 
of the cell here is W = 624. The variation in the scaled density 
fluctuation can be seen in Fig. 0. 
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FIG. 2: The steady state variation of the scaled density fluctuation 
N(x) across the width of the cell is shown here for two different 
cell widths. The parameters used for this simulation are: Na — 1, 
e = 0.75, and zy° = 0.077. A snapshot of the first plot is shown in 
Fig. ^ The dense regions coexist with dilute regions on both sides 
and are not necessarily restricted near the walls of the container, as 
shown in the second plot. 



distribution, whereas in the dense region the distribution func- 
tion shows non-Gaussian tails. Such non-Gaussian distribu- 
tions have been found in many other granular systems where 
other kinds of clustering phenomena have been found [ [l^ ] . 

It is to be noted that the two phase state is not observed if 
a constant temperature condition is used at the bottom wall 
for the same parameter values as those for which phase sepa- 
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FIG. 3: Pair distribution function in the dilute and dense phases. The 
static structure of the dilute phase is similar to a gas and the dense 
phase is similar to a liquid. 
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FIG. 4: The distribution function of horizontal velocities in the di- 
lute and the dense regions is plotted here. A Gaussian curve corre- 
sponding to the same local temperature T is also plotted along for 
comparison. It is seen that for the dilute region distribution function 
is very close to the Gaussian curve, whereas for the dense region the 
distribution function shows a non-Gaussian behavior. 



ration is observed with a vibrating wall. The physical reason 
for this is as follows. When a constant temperature condi- 
tion is enforced at the bottom wall, the mean kinetic energy 
of the particles after collision with the wall is independent of 
the kinetic energy before collision. This boundary condition 
tends to damp out temperature variations in the horizontal di- 
rection, since the post-collisional mean kinetic energy of the 
particles is independent of horizontal position. When a real 
vibrating wall is used, however, the collision of a particle with 
the wall only augments the kinetic energy of a particle, and 
so horizontal temperature variations are not damped out at the 



wall. This leads to the coexistence of regions with two differ- 
ent temperatures for a system driven by a vibrated surface for 
conditions under which a system with a constant temperature 
surface does not show this coexistence. 

The simulations also indicate that the horizontal position of 
the interface between the high and low density regions fluctu- 
ates in time, as shown in Fig. ||. To study the dynamics of the 
interface, we have plotted the Fourier transform of the density 
fluctuation at the interface. In the second plot of Fig. ||, we see 
that the dominant peak occurs only in the interface region and 
is clearly characterized by a single frequency, the correspond- 
ing time period (1.6 x 10^) matches the order of fluctuations 
seen in the time domain. 

It is natural to attempt to provide an explanation similar to 
[|^] for the coexistence of two densities in a vibrated bed sepa- 
rated into two compartments with a hole in the wall. However, 
it is shown here that such an explanation is not valid for the 
present case. In the spirit of [^, let us consider the possibility 
of a distribution of the particles into two "compartments" with 
two different densities (number of particles per unit width) 
Ni < N2. Coexistence requires that the horizontal fluxes 
of particles Fi and F2 between these two compartments be 
equal. The width Wi of the first compartment can be deter- 
mined from the constraint on the total number of particles. 
The constraints for a steady state distribution are 
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Fi = F2 
N2{W-Wi) = NW 



(2) 
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The above condition requires that the total horizontal flux 
F{N) be a non-monotonous function of N. This can be real- 
ized when the exchange between the two compartments takes 
place only at a particular height, as shown in [g] . When there 
is no wall separating the two compartments, the total flux of 
particles across any vertical surface is F{N) = N \^/27t. 
Since T - l/AT from Eq. (|l|), we have F - VN which is 
a monotonic function of N. Therefore, a low density the- 
ory cannot explain the existence of a steady distribution. The 
flux at high density has to be obtained using the density pro- 
files computed numerically by including the virial equation 
of state for the pressure — density relationship instead of the 
ideal gas equation of state [ p3| ] . It turns out that the variation 
of the fluxes for higher densities increases faster than a/TV. 
We have also computed the horizontal fluxes from simulation, 
by calculating the rate of particle collisions with the side walls 
for a system with different N and keeping all other parame- 
ters at a constant. For this simulation the width of the system 
is taken to be smaller than the length scale of separation so 
that the bed does not develop any inhomogeneity in the hori- 
zontal direction. The simulation results also indicate that the 
flux increases faster than a/TV in agreement with the theoreti- 
cal predictions. We have also verified that the simulations, in 
the case of the phase separation, do not violate the net zero 
flux condition across different vertical surfaces. This implies 
that the coexistence of two phases cannot be understood using 
simple steady state flux balance models such as the one pre- 
sented above, in contrast to the clusters (or phase separation) 
of [^, which could be explained using the dilute bed kinetic 
theory of [R] . However, it can be shown [pj]] that the variation 
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of the density and the temperature in the homogeneous state 
just before the occurrence of this instabihty can be adequately 
represented by the kinetic theory of [|l4|], where the effects of 
the periodicity of wall vibration are neglected. This indicates 
that the stability of the system can be understood in terms of 
simple equations. 



50 




(D xlO" 



FIG. 5: Dynamics of the interface. A close inspection of the local 
scaled density fluctuations N{x) in top plot shows that the cluster 
interface is not static but varying with time. The second plot shows 
the Fourier transform (FFT) of the density fluctuation in the time 
domain at different horizontal positions. A peak is found at x = 
770 near the interface region. The characteristic radian frequency is 
around 4 x 10~^, which corresponds to a time period of 1.6 x 10^ 
roughly of the order of interface fluctuation in the top plot. 



To conclude, a spontaneous symmetry breaking phe- 
nomenon is observed in a vibrated bed in a regime where it 
has not been observed before, i.e., in the dilute limit where 
the area fraction is small, the amplitude of the vibrating sur- 
face is small compared to the mean free path and the fre- 
quency is large compared to the frequency of mean motion 
of the bed. This appears similar to that reported by [||], but a 
detailed analysis indicated that the present phenomenon can- 
not be explained the simple model used there. This requires 
a more complicated model where the fluctuating nature of the 
interface and possibly the non-Gaussian nature of the distri- 
bution function in the dense phase is incorporated. A stability 
analysis of this state could have implications in understanding 
instabilities in driven granular materials and in pattern forma- 
tion in vibrated beds. 
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